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Abstract 



This paper is concerned with adaptive signal control problems on a road network, 
using a link-based kinematic wave model (Han et al. , 2012). Such a model employs the 



Lighthill-Whitham-Richards model with a triangular fundamental diagram. A variational 



type argument (Laxl 1957| Newell, 1993) is applied so that the system dynamics can be 



determined without knowledge of the traffic state in the interior of each link. A Riemann 
problem for the signalized junction is explicitly solved; and an optimization problem is 
formulated in continuous-time with the aid of binary variables. A time-discretization 
turns the optimization problem into a mixed integer linear program (MILP). Unlike the 
cell-based approaches (Daganzo 1 1995 Lin and Wang, 2004 Lo, |1999b| ), the proposed 



framework does not require modeling or computation within a link, thus reducing the 
number of (binary) variables and computational effort. 

The proposed model is free of vehicle-holding problems, and captures important fea- 
tures of signalized networks such as physical queue, spill back, vehicle turning, time- varying 
flow patterns and dynamic signal timing plans. The MILP can be efficiently solved with 
standard optimization software. 



1 Introduction 



Traffic signal is an essential element to the management of the transportation network. For 
the past several decades, signal control strategies have evolved from ones developed based on 
historical information, often referred to as the fixed timing plan, to the generation of control 
strategies in which the control system is fully responsive. In the latter case, the cycle lengths 
and splits of the signal are determined based on real-time information. Representatives of 



2000), SCAT (Sims and Dobinson, 1980) and SCOOT (Hunt et al., 1982). 



such signal-control systems are OPAC (Gartner, 1983), RHODES (|Mirchandani and Head 



The performance of a traffic signal control system depends on the optimization proce- 
dure embedded therein. We distinguish between two optimization procedures: 1) heuristic 
approach, such as those developed with feedback control, genetic algorithms and fuzzy logic; 
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and 2) exact approach, such as those arising from mathematical control theory and mathe- 
matical programming. Among these exact approaches, the mixed integer programs (MIPs) are 



of particular interest and has been used extensively in the signal control literature. |Improta 



and Cantarella (1984) formulated the traffic signal control problem for a single road junction 



as a mixed binary integer program. Lo (1999b) and Lo (1999c) employed the cell transmission 



model (CTM) ( |Daganzo , 1994, 1995) and casted a signal control problem as mixed integer 
linear program. In these papers, the author addressed time-varying traffic patterns and dy- 



namic timing plan. In Lin and Wang (2004), the same formulation based on CTM was applied 



to capture more realistic features of signalized junctions such as the total number of vehicle 
stops and signal preemption in the presence of emergency vehicles. One subtle issue asso- 
ciated with CTM-based mathematical programs is the phenomena known as traffic holding, 
which stem from the linear relaxation of the nonlinear dynamic. Such an action induces the 
unintended holding of vehicles, i.e., a vehicle is held at a cell even though there is capacity 
available downstream for the vehicle to advance. The traffic holding can be avoided by in- 
troducing additional binary variables, see Lo (1999a). However, this approach ends up with 



a substantial amount of binary variables and yields the program computationally demanding. 
An alternative way to treat holding problem is to manipulate the objective function such that 
the optimization mechanism enforces the full utilization of available capacities in the network. 
This approach however, strongly depends on specific structure of the problem and the under- 



lying optimization procedure. Specific discussion on traffic holding can be found in Shen et 



al. (2007) 



This paper presents a novel MILP formulation for signal control problem based on the 



link-based kinematic wave model (LKWM). This model was proposed in |Han et al. (2012) as 
a continuous-time extension of the Lighthill-Whitham- Richards (LWR) model (Lighthill and 



Whitham, 1955 Richards, 1956) to networks. This model describes network dynamics with 



variables associated to the entrance and exit of each link. It employs a Newell-type variational 



argument (Newell, 1993; Daganzo, 2005) to capture shock waves and vehicle spillback. Ana- 



lytical properties of this model pertaining to solution existence, uniqueness and well-posedness 



are provided in Han et al. (2012). A discrete-time version of the LKWM, known as the link 



transmission model, was discussed in Yperman et al. (2005). In contrast to the cell-based math 



programming approaches where the variables of interest correspond to each cell and each time 
interval, the model proposed in this paper is link-based, i.e. the variables are associated with 
each link and each time interval. The resulting MILP thus substantially reduces the number 
of (binary) variables and hence the computational effort. In addition, the link-based approach 
prevents vehicle holding within a link without using binary variables. The resulting mixed 
integer linear program can be solved efficiently with commercial optimizers such as CPLEX. 

The formulation in this paper captures key phenomena of vehicular flow at junctions 
such as the formation, propagation and dissipation of physical queues, spill back and vehicle 
turning. It also considers important features of signal control such as dynamic timing plan and 
time-varying flow patterns. In the remainder of this introduction, we briefly review the LWR 
model and the variational method for solving the corresponding Hamilton-Jacobi equation. 
This will serve as preliminary background as we proceed in Section [2] to discuss our link-based 
network mode. 



1.1 Lighthill-Whitham- Richards model 



Following the classical model introduced by Lighthill and Whitham (1955) and Richards 



(1956), we model the traffic dynamics on a link with the following first order partial dif- 
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ferential equation (PDE), which describes the spatial-temporal evolution of density and flow 

(1.1) 



9 p(t, x) + -^-f(p(t, x)) 



dt 



where p(t, x) : [0, +oc) x [a, b] — » [0, pj] is average vehicle density, f(p) : [0, fP am ] [0, C] 
is average flow. fP am is jam density, C is flow capacity. The function /(•) articulates a 
density-flow relation and is commonly referred to as the fundamental diagram. 



Classical mathematical results on the first-order hyperbolic equations of the form (1.1) can 



be found in Bressan (2000). For a detailed discussion of numerical schemes for conservation 
laws, we refer the reader to Godunov (1959) and LeVeque (1992). A well-known discrete 



version of the LWR model, the Cell Transmission Model (CTM), was introduced by Daganzo 



(1994, 1995). PDE-based models have been studied extensively also in the context of vehicular 



networks, with a list of selected references including Bretti et al. (2006); Coclite et al. (2005) 



(2003); Lebacque and Khoshyaran (1999, 2002) 



Daganzo (1995); Herty and Klar (2003); Holden and Risebro (1995); Jin (2010); Jin and Zhang 



1.2 The Hamilton- Jacobi equation and the Lax-Hopf formula 



Initially introduced by Lax (1957, 1973), then extended in Aubin et al. (2008) and LeFloch 



(1988), and applied to traffic theory in Claudel and Bayen (2010); Daganzo (2005); Friesz et 



al. (2012); Han et al. (2012), the Lax-Hopf formula provides a new characterization of the 



solution to the hyperbolic conservation law and Hamilton-Jacobi equation. The Lax formula 
is derived from the characteristics equations associated with the H-J equation, which arises 
in the classical calculus of variations and mathematical mechanics. The reader is referred to 



Evans (2010) for a detailed discussion. 



Let us introduce the function 7V(-, •) : [0, +oc) x [a, b] —> R, such that 



d_ 
dx 



N{t, x) 



p{t, x) 



d_ 
dt 



N(t, x) = /(p(t, x)) 



(1.2) 



The function N(t, x) is sometimes referred to as Moskowitz function or Newell-curves. It has 
been studied extensively, in Claudel and Bayen ( |2010 ); Daganzo (2005); Moskowitz (1965); 
Newell (1993). A well-known property of N(-, •) is that it satisfies the following Hamilton- 
Jacobi equation 



d_ 
di 



N(t, x)-f 



d_ 

dx 



N(t, x] 



(t, x) e [0, +oo) x [a, b] 



(1.3) 



[a, 6], and c(-, •) : Q —> 



Let ^ be a subset of [0, +00 

the value of 7V(-, •) on £1. The solution to the H-J equation (|1.3[) with condition c( 



be a value condition which prescribes 

is given 



by the Lax-Hopf formula |Claudel and Bayen| ( |20To| ) ; |Daganzo| ( |2005| ) ; |LeFloch| ( fl988| ) 



N(t,x) = m{{c(t-T,x-Tu)+T f(u)} (1.4) 

such that u e [f{p jam - ), /'(0+)], T > and (t - T, x - Tu) e ft. Where /*(•) is the 
Legendre (concave) transformation of /(•). 

1.3 Oganization 

The rest of this article is organized as follow. In section [2j we present a network model 
known as the link-based kinematic wave model ( Han et al| |2012 ). Section [3] formulates the 
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traffic signal control problem in both continuous- and discrete-time, based on the LKWM. 
The discrete-time problem is further formulated as a mixed integer linear program. Section [I] 
presents a numerical example, which demonstrates and evaluates the proposed formulation. 



Link-based Kinematic Wave Model 



In this section, we present a kinematic wave model on networks, with a triangular fundamental 
diagram for each link. Unlike the cell-based models |T)aganzo (1994, 1995), the proposed model 
does not require modeling or computation in the interior of the link. For this reason, we call 
it the link-based kinematic wave model. 



2.1 State variables of the system 

Consider the link represented by an interval [a, 6], with b — a = L > 0. In the derivation 
of the LKWM, we select flow x) and regime r(t, x) as the state variables for the link, 
instead of density. It is obvious that a single value of q corresponds to two traffic states: 
1) the free flow phase (r = 0); and 2) the congested phase (r = 1). Therefore, the pair 
x), r(£, x)) E [0, C] x {0, 1} determines a unique density value. This simple observation 
gives rise to the following map 

^(■) : [0, C] x {0, 1} -+ [0, (f am \ (q, r) ^ p (2.5) 



2.2 Riemann problem at a junction with one incoming link 



Extension of the kinematic wave model to a network turns out to be subtle; the issues asso- 
ciated therein are 1) a proper definition of a weak entropy solution at a junction of arbitrary 
topology; 2) uniqueness and well-posedness of the entropy solution. The reader is referred to 
Han et al.| ( |2012[ ); |Garavello and Piccoli| Q2006P ; [jm| fl2010P ; |Jin and Zhang] fl2003| ) for some 
specific discussion. A junction model can be analyzed by considering a Riemann problem, 
which is an initial value problem with constant datum on each incoming and outgoing link. 
Due to space limitation, instead of a comprehensive discussion on various types of Riemann 
problems, we focus on the Riemann problem for a particular junction, that is, the one with 
one incoming link and n > 2 outgoing links. This is because we assume that during one signal 
phase, cars from only one incoming link can enter the junction. 

In order to model vehicle turning, we fix a traffic distribution matrix 

A = (a lj2 «i 5 3 . . . ai, n +i) 

where < ot\^ < 1, i = 2, . . . , n + 1, Y^l=2 a M ~ 1- ^ ne coefficients ol\^ determines how 
the traffic from the incoming link 1\ distributes in percentages to the outgoing link 1{. For 
simplicity, we assume A is time-independent. Note that there is no substantial difficulty with 
transforming our modeling framework to deal with time-varying distribution matrices; such 
extension, however, requires additional information on route choices, which is beyond the 
scope of this paper. 

The next theorem characterizes the solution to the Riemann problem at junction with one 
incoming link. 

Theorem 2.1. Consider a junction with one incoming link I\ and n > 2 outgoing links 
I2, . . . , I n +i- For every initial data yi,o 5 • • • 5 2/n+i,o £ [0, C] x {0, 1} ; there exists a unique 
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n + 1 -tuple 

VI, £n+l G [0, C] X {0, 1} 

where yi — (ft, fi), such that the solutions to the initial-boundary value problems at the junction 



m p{t,x) + ^h(p{t,x)) = 



p(0, x) 



, § i p(t,x) + ^f i (p(t,x)) = o 
p(0, x) = ^(^, ) 



n+1 



Z5 t/ie admissible weak solution to the junction problem in the sense defined in \Coclite et al. 



(2005). In addition, we have the following characterization: the boundary states (ft, f^), i 
1, . . . , n + 1 are given &y 



91 



mm < g^ aa \ 



„max 
92 



„max 
93 



ri 




«1,2 «1,3 

if rift 
if rift 



9n+l 



1 

0, 9i = 9i,o 
0, 9i < 9i,o 
n+1 





if r ifi = 1, ft = 9i,o 
i/ H,o = 1,9* < 9i,o 



where 



9;,o + ^,o(C; -9i,o), i = 1 

Ci + riflfaifl -Ci), i = 2, 



n + 1 



(2.6) 



(2.7) 



(2i 



(2.9) 



(2.10) 



Remark 2.2. T/ie quantity qV 1 ^ is the maximum flux an incoming (outgoing) link can send 
(receive) - a quantity identified as demand (supply) by Lebacque and Khoshyaran\ $199£\ 2002). 



The aforementioned map (yi } o)i=i,...,n+l ^ (&)*=!.,.. .,n+l is commonly referred to as the 
Riemann solver, see |Garavello and Piccoli (2006) for a formal discussion. Theorem 2.1 de- 
scribes the Riemann solver using the new state variables q and r. The verification of theorem 



2.1 is straightforward. For junctions with arbitrary topology, the Riemann Solvers are not 
available in closed- form. Yet, in the case of one incoming link, we are able to express the 
Riemann solver explicitly. 

In the next subsection, we analyze shock formation and propagation within one single link. 
The location of the shock wave is crucial as it determines the regime variable r associated with 
the two boundaries of the link. 



2.3 Shock formation and propagation within the link 

We focus on solutions generated by assuming an initially empty network, i.e. y^o — (0, 0). 
The key to our analysis is the location of a so-called separating shock, which divides each link 
into two zones: free flow zone (r = 0), and congested zone (r = 1). We begin with the fact 
that if the network is initially empty, then there can be at most one separating shock on each 
link 
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Lemma 2.3. For every link I{ and any solution yiit, x) — {qi{t, x), r^(t, x)) with yi(0, x) — 
(0, 0), the following statement holds: 

1. For every t > 0, there exists at most one x*(t) E (ai, bi) such that r^(t, x*(t)—) < 

n(t, x*(t)+) 



2. For all x E [a^, 



ri(t, a) = 0, if x < x*(t) 
ri(t, x) = 1, if x > x*(t) 



Proof See Bretti et al. (2006). 



□ 

According to Lemma 2.3, the separating shock emerges from the downstream boundary 
b{ of the link, and propagates towards the interior of the link. The speed of this separating 



shock is given by the Rankine-Hugoniot condition Evans (2010). It is clear that as long as the 



separating shock remains in the interior of the link Ii, the upstream and downstream boundary 
conditions do not interact. Thus the exit of a link remains in the congested phase; while the 
entrance remains in the free flow phase. Consequently, the Riemann Solver in Theorem 2.1 is 
expressed entirely with exogenous parameters C\ and ol\^. On the other hand, if the separating 
shock reaches either boundary, it becomes a latent shock. Two cases may arise. 

i) The the shock reaches the exit, i.e. x*(t) = b{. In this case, the current link is dominated by 
free flow phase. The boundary condition at x = ai directly influences the boundary condition 
at x = bi, in a way expressed by 



Qi (*, bi-) 



Li 

t- — , ai 



(2.11) 



where Li is the length of the link, k\ is the speed of forward wave propagation. See Figure [T] 
for an illustration. 

ii) The shock reaches the entrance, i.e. x*(t) = a^. In this case, the current link is dominated 
by congested phase. The boundary condition at x = bi directly affects the boundary condition 
at x — ai 

Li 



qi (t, di+) = qi yt 
where W{ is the speed of backward wave propagation. 



Wi 



(2.12) 



In either case, the Riemann solver involves boundary flows (2.11), (2.12), which are en- 



dogenous. The next key step towards the link-based flow model is the detection of these two 
extreme cases. This can be done with the Lax-Hopf formula. 



2.4 The variational approach for detecting latent shock 

This section provides sufficient and necessary condition for the occurrence of the latent shock. 



The derivation is omitted for brevity, we refer the reader to Han et al. (2012) for a detailed 



discussion. Define for each link Ii the cumulative entering and exiting vehicle numbers 



i(s, a^ ds, 



Ni 



i.down 



qi(s, bi)ds 



io Jo 

Recall that the separating shock x*(-) : [0, +oc) —> [ai, bi] is a continuous curve in the t — x 
domain. The following theorem provides sufficient and necessary condition for the occurrence 
of the latent shock. 
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case i) case ii) 

Figure 1: Example of latent separating shocks. Case i): the separating shock reaches the right 
(downstream) boundary. Case ii): the separating shock reaches the left (upstream) boundary. 



Theorem 2.4. Let A^(-, •) : [0, +oc) x [a^, bi] be the unique viscosity solution to the Hamilton- 
Jacobi equation {l.Sfy , satisfying zero initial condition, upstream boundary condition Ni^ up (-) 
and downstream boundary condition Ni^own{')- Then for all t > 0, 



x^ (^) — 
x*(t) = bi 



Ni, up CO > N hdown [ t - ^ ) + pf m L % 



Ni,up [ t 

l%i 



^ Ni^down (^) 



(2.13) 
(2.14) 



Remark 2.5. The significance of criteria (2.13)-(2.14) is that the two extreme cases can 



be detected without any computation within the link. This is because N up (-), N^owni') are 



determined completely by the boundary flows. Theorem \2.J\ is the key ingredient of the LKWM, 
which allows the network model to be solved at the link level. 

Analytical properties of the LKWM pertaining to solution existence, uniqueness and well- 
posedness are provided in Han et al. {2012 ). 



3 Traffic signal control problem based on the LKWM 

In this section, the signal control problem is formulated with LKWM. We start with a single 
junction, with two incoming links ii, I2, and two outgoing links ^3, I4 (Figure [2]). Each link 



1 



2 



Figure 2: A signalized junction with two incoming links and two outgoing links, 
is represented by a spatial interval [a^, i — 1, 2, 3, 4. The fundamental diagrams for each 



7 



link is given by 



Hp) - \ kip Pe[0,p*] z -i234 

[-Wi{p-pi ) P^KPnPi \ 

with d = ki p\ being the flow capacity. We make the following assumptions 
Al The network is initially empty. 

A2 Drivers arriving at the junction distribute on the outgoing roads according to some known 
coefficients: 

/ «l,3 OL\A 



A 

where denotes the percentage of traffic coming from link li that distributes to outgoing 



link Ij. 



In the problem setting, the flows %(•), i = 1, 2 entering links I\ and I2 are known. In 
practice, qi{-) can be measured at the entrance using fixed sensors such as loop-detectors. 



3.1 Continuous-time formulation 

In this section, we formulate the constraints of the system in continuous time. In Section [3. 2| 
we will reformulate the system dynamic as linear constraints in discrete time using binary 
variables. Notice that the discussion in this section can be the building block for extension to 
networks with multiple intersections. 

Let us fix the planning horizon [0, T] for some fixed T > 0. Introducing the piecewise- 
constant control variables i^(-) : [0, T] —> {0, 1}, i = 1, 2, with the agreement that Ui(t) = 
if the light is red for link 1^ and U{(t) = 1 if the light is green for link 1{. It is convenient to 
use the following set of notations. For i = 1, 2, 3, 4. 





the 


flow of cars entering link 7^, 




the 


flow of cars exiting link 7^, 




the 


binary variable that indicates the regime at x = 




the 


binary variable that indicates the regime at x = b{ — 


-max ( m \ 
Hi \ )•> 


the 


maximum flow allowed to enter the link Xj, 


"i \ ' ) •> 


the 


maximum flow allowed to exit the link Zj, 




the 


cumulative number of cars that have entered link ^, 


Ndown,i{') i 


the 


cumulative number of cars that have exited link 7^, 


V>i(')i 


the 


signal control variable for link Zj, 



Theorem 3.1. The dynamics at the junction (FigureWfi with signal control can be described 
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by the following system of differential algebraic equations (DAE) with binary variables. 

^N up>i (t) = qi(t), ^AWm,i(*) = q_i(t), 

'l, if N up ^t)>N down ^ i {t-^)+fy i ^ m L i 



1, 2, 3, 4 



n{t) 
hit) 



0, otherwise 



0, if N uPil yo - 

1, otherwise 



^ Ndown,i{f) 



-max ( 



(t) = C i + f i (t)[q i (t-^)-Q 



Wi 



iax (t) 



0, if Ui (t) = 

~ rx{t)i } ,/ Ui(i) = i 

q k (t) = ai tk qi(t) + a,2,kfa{t), 



1, 2, 3, 4 

1, 2, 3, 4 

1, 2 
3,4 

1, 2 



fc = 3, 4, 



ui(t) + u 2 (t) = 1 for all t € [0, T] 



(3.15) 
(3.16) 

(3.17) 

(3.18) 
(3.19) 

(3.20) 

(3.21) 
(3.22) 



Proof. (3.15) is by definition. For % = 1,2,3,4, if the separating shock on fink ij reaches the 



entrance (exit b{), then the regime variabie fj = 1 (fj = 0). Then (3.16)-(3.17) foilows from 
Theorem 12.41 

The demand function qfi iax {-) for incoming links and the supply function <^" aa: (-) for out- 



going links are given by (2.10): for i = 1, 2, if fi(t) = 1, then g[™ ax (t) = Q; otherwise if 



ri(t) — 0, then according to (2.11), 

! (t) = 6i 



~max ( 



U 

t- — , a, 
Ki 



hi 



This shows (3.19). One can similarly show (3.18) using (2.10) and (2.12). 



For z = 1, 2, if 'Ui(t) = which means the light is red, then the flow allowed through is 



zero, otherwise, it is given by(2.6). This proves (3.20) 



(3.21) follows from the definition of the splitting parameters a^, i = 1, 2, k = 3, 4. (3.21) 
guarantees that at each time, there are one and only one incoming road that has green light. 

□ 



3.2 Discrete- time formulation 

In this section, we present the discrete-time version of the optimization problem in Theorem 



3.1| Let us introduce a few more notations for the convenience of our presentation. Consider 
a uniform time grid 

= t° < t 1 . . . < t N = T, t j - t j ~ l = St, j = 1, . . . , N 

Throughout the rest of this article, we use superscript c j' to denote the discrete value evaluated 
at time step t j . In addition, we let U/ki = A{St, U/wi = A\8t, A{ e N, A\ e N, i = 
1,2,3,4. 

Approximating the numerical integration with rectangular quadratures, we write equality 



(3.16) and (3.17) in discrete time as 
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st J2 i-stYti+^u < M(i-ff) 



< 



fc-Aj 

stj: 

3=0 



A? < k < N, i 



1, 2, 3, 4 



StYi + pi^Li > -Mff + e 

3=0 



1 % 



3=0 
k-A{ 



3=0 
k 



(3.23) 

A{ < k < N, i = 1, 2, 3, 4 (3.24) 



3=0 3=0 

where ff^rf G {0, 1}. .M G is a sufficiently large number, £ E is a sufficiently small 
number. Constraints ( |3.23 ) and (3.24) determines the regime variables associated with the 
two boundaries of each link. Once the flow phases are determined, the demand and supply 



functions (3.18), (3.19) are re- written in discrete time as 



M(l-fj) < q™ ax ' 3 < C i +M(l-r, 



d + M (rj - 1) < C ' J < C 



s 



-Mr 3 - < q, 



max, j 



< 



f 



+ Mr 3 , 



1, 2, 3, 4 (3.25) 



i = 1, 2, 3, 4 (3.26) 



Next, let us re-formulate (3.20). Introducing dummy variables Ci ? C2 ' — ^ — ^ ' such that 

-max,j -max, j ^ 



Ci = mm <^ q 



max, j 



chlz aid 



i = 1, 2 



(3.27) 



Then the discrete-time version of (3.20) can be readily written as 
f0 < ^ < A^^ 

l^+^K-i) < g? < CJ ' ' J 



(3.28) 



In order to write ( 3.27[ ) as linear constraints, one could write it as three "less or equal" 
statements, which is simple but bear the potential limitation of traffic holding. Instead, one 
may introduce additional binary variables ^, r]j and real variables f3j for i = 1, 2, j = 1, . . . , JV, 



such that (3.27) can be accurately formulated as 



-max,i I k a i-3 ^ qj ^ -max,j / 

-max,j I k a f i ti\ ^ oj ^ -max A / 

Pi-M{i-4) < ci < p{ 

Finally, we have the obvious relations 

q{{t) = ai,Jfe9i(*) + a 2,kQ%(t) k = 3, 4, 



i = 1, 2 



(3.29) 



1 N 



(3.30) 
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and 



u\ + uL — 1 



3 = 1,...,N 



(3.31) 



The proposed MILP formulation of signal control problem is summarized by (3.23 )-( 3.26| ) 
and (3.28)-(3.31). This formulation captures many desirable features of vehicular flow on 



networks such as physical queues, spill back, vehicle turning, and shock formation and prop- 
agation (although not explicitly). The signal control allows time- varying cycle length and 
splits, as well as the utilization of real-time information of traffic flows. 



3.3 Bound the separating shock 

At the end of this section, we discuss an additional linear constraint that ensures that the 
congested phase on each link is bounded. Such condition is closely related to travel delay: if 
the congested phase remain bounded, there will be a reasonable upper bound for the travel 
time of each driver. Let us consider a single link [a, b]. If we wish to bound the separating 
shock within the interval [c, b] for some a < c < 6, this means (with the same notation as 
before) that 

rt r t-^ 

/ q(s, c)ds < W q(s) ds + fP am (b - c) (3.32) 

Jo Jo 

(3.32) follows by applying ( 2.13| ) to the interval [c, b\. It is helpful to notice that since [a, c 
remains in the free flow phase, q(s, c) must be equal to q (s — ^jr) 
the congested phase to remain in [c, b] becomes 



Thus the condition for 



< N down [t-—)+p> am {b-c) 



w 



(3.33) 



Condition (3.33) can be easily written as linear constraint in discrete time. If one wishes to 



include (3.33) in the objective function instead of using it as a constraint, he/she may simply 



minimize the difference between the left and right hand sides of (3.33) 



4 Numerical example 

In this section, we consider the simple network consisting of two junctions, as shown in Figure 
[3j We will show the numerical result of optimal signal control at this junction obtained by 
the MILP summarized in the previous section. 



1 v 



2v 



6v 



7 v 



Figure 3: Test network with two signalized intersections. 
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4.1 Numerical setting 

We assume that the fundamental diagrams for all seven links are the same, that is, for i = 
1,...,7 

k{ = 30 mile/hour, W{ = 10 mile/hour, p 3 ^ 171 = 400 vehicle/mile, C{ — 3000 vehicle/hour 

The lengths of all links are set equally to be 0.3 miles. We choose a time grid of 20 intervals 
and a time step of 0.005 hour (18 seconds). The flow entering links I\,l2,H are chosen to 
be time-varying functions whose value at each time interval is randomly generated between 
veh/hour and 3000 veh/hour. In addition, to ensure the performance of signal control, we 
include further constraint that the congested phase must not pass the mid-point of the link, 
i.e. it remains on the spatial interval [0.15,0.3]. This is done by invoking constraint (3.33). 



The MILP is solved with ILOG Cplex 12.1.0, which runs with Intel Xeon X5675 Six-Core 
3.06 GHz processor provided by the Penn State Research Computing and Cyberinfrastructure. 



4.2 Numerical results 

The solution time of the MILP described above is 0.49 seconds. In order to have a clear 
visualization of the optimal signal strategy and the separating shocks on each link, we use 
the boundary datum q^ qi, i = 1, 2, 3, 4 obtained from the optimal solution to construct 



solutions to the Hamilton- Jacobi equation (1.3) using Lax-Hopf formula. For the H-J equation, 
the separating shock no longer represents discontinuity, rather, it is displayed as a 'kink' 
(discontinuity in the first derivative). The Moskowitz functions N(t, x) for links ^3, I4 are 
shown in Figure [4] and [5j respectively. The Moskowitz functions for link I\ and I2 viewed from 
a different angle are shown in Figure [6] and [7| 

One can clearly observe, in each figure, two types of characteristics lines: forward ones 
and backward ones, representing the free flow and congested phase, respectively. The shared 
boundary of the two regimes is precisely the separating shock wave, which we managed to 
implicitly handle with variational method. It is also clear that the congested region never 



crosses the middle point of the link throughout the planning horizon, thanks to (3.33). 




t (hour) 




x (mile) 



t (hour) 



x (mile) 



Figure 4: Moskowitz function for link 1%. Figure 5: Moskowitz function for link I4. 



12 



0.15 

x (mile) 



0.15 

x (mile) 



Figure 6: Moskowitz function for link I\. Figure 7: Moskowitz function for link I2. 

5 Conclusion 

This paper proposes a link-based signal control problem. The key ingredient of our formulation 
is the use of the Lax-Hopf formula to detect latent shocks and hence to identity the regimes 
(free flow/congested) at the entrance and exit of each link. The analytical framework allows 
us to further bound the congested region and to avoid spillback. The problem of adaptive 
signal control is formulated in discrete time as a mixed integer linear program. The resulting 
MILP requires fewer (binary) variables compared to cell-based approaches. 

One limitation of the current approach is the lack of computational tractability, when 
the problem size scales up. Future research aims at exploring heuristic optimization algo- 
rithms such as the genetic algorithms; and new computational paradigms involving parallel 
computing. 
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